Full Paper 



AMAE Int. J. on Manufacturing and Material Science, Vol. 3, No. 1, May 2013 



Effect of Piezoelectric Layer on Beam Parameters 

using Zigzag Theory 

Mohammad Sarim 1 , Dr. M. Naushad Alam 2 
1 Post Graduate Student, Dept. of Mechanical Engineering, 
ZHCET, AMU, Aligarh, India 
Email: mohammadsarim@zhcet.ac.in 
2 Associate Professor, Dept. of Mechanical Engineering, 
ZHCET, AMU, Aligarh, India 
Email: naushad7863@rediffmail.com 



Abstract — An efficient higher order theory is presented for 
static analysis of multilayered composite beams with 
piezoelectric layers embedded or bonded to the surface, under 
static electromechanical load. In this theory, the in-plane 
displacement field is taken as a combination of a layer-wise 
linear variation and a cubic variation across the thickness. 
Transverse normal strains are neglected. The electric field is 
also approximated as piecewise linear across the sub layers. 
The displacement field is expressed in terms of only three 
primary displacement variables excluding electric potential 
variables by enforcing the conditions of zero transverse shear 
stress at top and bottom of the beam and its continuity at layer 
interfaces under general electromechanical loading. The effect 
of thickness of the piezo-layer is observed for various loading 
conditions. Also, the effect of beam lay-up on various 
parameters is studied. 

Index Terms — zigzag theory, hybrid beam, piezoelectric layer, 
static analysis, thickness, electromechanical loading. 

I. Introduction 

Hybrid beams and plates are made of an elastic substrate 
with piezoelectric layers embedded or bonded to its surface. 
Due to inhomogeneity in mechanical properties across the 
thickness and presence of piezo-layers, analysis of these 
structures requires appropriate electromechanical modeling. 
3D analytical solutions have been presented for the piezo- 
electric response for simply supported rectangular plates in 
Heyliger, 1994 [1]; Lee and Jiang, 1996 [2]. Early works on 2D 
theories for hybrid laminated beams were based on classic 
laminate theory approximation for mechanical field without 
considering the coupling between mechanical and electric 
fields (Tzou, 1989 [3]). The limitation of CLT not including 
shear deformation effect was addressed by First Order Shear 
Deformation Theory (FSDT) (Kapuria et al., 1997b [4]) and 
refined third order theory of Reddy (1984) [5]. 

This paper presents an efficient higher order zigzag theory 
for the static analysis of hybrid piezoelectric beams subjected 
to electromechanical loads. The displacement field is taken 
as a combination of a linear variation layerwise and a cubic 
variation across the thickness. Transverse normal strains are 
neglected. The displacement field is expressed in terms of 
only three primary displacement variables excluding electric 
potential variables by enforcing the conditions of zero trans- 
verse shear stress at top and bottom of the beam and its 
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continuity at layer interfaces under general electromechani- 
cal loading. The governing equations and boundary condi- 
tions are derived from the principle of virtual work and ana- 
lytical Navier solution is obtained using Fourier series expan- 
sion. The effect of thickness of piezoelectric layer is examined 
for various beams keeping the total thickness of beam same in 
all cases and varying the thickness of piezo-layer. The results 
are plotted for various beam span-to-thickness ratios. The 
effect of beam lay-up on various beam parameters is also stud- 
ied and the results are given for cross-ply and angle -ply 
beams. 

II. Formulation Of Theory 

A. Geometry of the Hybrid Beam 

Consider a hybrid beam of solid cross-section of width b 
and thickness h, consisting of orthotropic layers with their 
principal material axis along the fibers at an arbitrary angle to 
the Cartesian coordinate axis-x along its length. The geometry 
of the beam is shown in Fig. 1 . Some of the layers can be 
piezoelectric with class mm2 symmetry and poling along the 
thickness axis-z. The origin of the coordinate system (x,z) is 
taken at the middle surface of the beam. The integer k denotes 
the layer number which starts from the bottom of the laminate. 
The distance from the reference plane to the bottom surface 
of the k lb layer is denoted by z . Each layer is perfectly 
bonded to its adjacent plies. 

B. Constitutive Equations 

For a piezoelectric continuum which exhibits class mm2 
symmetry with respect to principal material axes x , x, jc,(= z) 
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Figure 1. Geometry of the hybrid beam 
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and is polarized along direction z, the constitutive equations 
referring to the orthogonal coordinate system (x,y,z) are given 
by (Auld, 1973) [6]; Tzou and Bao, 1995 [7] as 

e = S a + d T E , })= J<j + ~E (1) 
where the superscript T denotes matrix transpose. S~ , d~ > 
and £ are respectively the transformed compliance 
coefficients, piezoelectric strain constants and constant stress 
dielectric constants. 

For a beam with small width, the following assumptions 
are made: 



o ,cr ,T ,T ,E =0 

z 7 y 7 yz 7 xy 7 y 
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On use of (2), the constitutive relations in (1) reduce to 
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C. Displacements and Potential Field 

The presented ID beam model combines a third order 
zigzag approximation for the in plane displacement u through 
the thickness of the laminate with a layerwise approximation 

for the electric potential <fi such that the continuity of the 
transverse shear stress at the interfaces between adjacent 
layers is satisfied in presence of the electromechanical field. 
The transverse displacement field w is obtained by 
superimposing a constant field across the thickness with a 
layerwise contribution that accounts for the piezoelectric out- 
of-plane normal strain induced by the electric potential. 
Accordingly, the displacement and the potential fields are 
assumed in the following form: 

u(x, z) = u k (x) - zw (x) x + zy k (x) 

—k N 

+Z z ^{x)+zrj{x)+a £/'(zty£ (4a) 



w(x, Z )=w {x)-d 33 jy j (z)<t> j 

H 

N 

</>(x,z)=^ x ¥ i (z)(l> i (x) (4b) 

with a =d 3 3-e l5 1 q" 55 and /'(«) = J o V ; (z)& .Here a sub- 
script comma denotes differentiation, u. and "F t denote dis- 
placement and rotation variables of the kth layer. N is the 
number of points i across the laminate thickness for describ- 
ing the variation of the potential field in the thickness direc 



tion. $ 1 are the electric potentials at points z' and >p * ( z ) are 
linear interpolation functions. 

The functions u p and v\ are determined using the 

following conditions: 

1. The top and bottom surfaces of the beam are traction-free 

1. e. r_ x |. -JW2 = . 

2. The shear stress is continuous across each layer interface. 

3. The displacement u is continuous across each layer 
interface. 

The final expression of u is obtained as: 



(5) 



where 



= < + zR[ + z 2 R, + z'R,, F t '(z) = a f i (z) + R^ 
It can be seen that above displacement field contains 
three primary variables (u Q ,w Q , >F , ) excluding the contribution 

from the potential variables ^ ' . The number of the primary 
variables is thus the same as in the FSDT 

D. Governing Equations 

The equilibrium equations and the variationally 
consistent boundary conditions can be formulated in a weak 
form using the principle of virtual work given by: 



\{u i 8s i +Z)<^ )dV =\(p i Su i -qS(/>)dA 



(6) 



where V and A represent the volume and the surface of the 
piezoelectric continuum. 

| v o.5s.dV and j v D.5(f).dV represent the virtual work 
done respectively by the internal forces and the electric 

field. J" A p.Su.dA and j" A qd(f>dA are the virtual external works 
done by the applied surface tractions p. and charge density q 
respectively, a . ,e. , u. and D denote respectively the 
components of stress, strain, displacement and electric 
displacement. 

The electromechanical equilibrium equation in terms of 
primary field variables is derived in Kapuria (2001) [8]. 

[L ij ]U = P (7) 



E. Static Analysis 

The boundary conditions at simply supported ends are: 
N x = w = M x = P x = s; =0, 

j= 1, ... ,N atx = 0,a. (8) 

The solution of the governing equation (7) which 
identically satisfies the boundary conditions as given by (8) 
is expressed in terms of Fourier Series as: 

(w, , (j> ' , N , , M x , P , Sj , G ' ) = 
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, $ , N x , M x , P , S J x , G' ) n sin nx 

n-l 

(« , % , Q x , Q x ,H') = . ¥i . Q , Q[ ,H')„cosnx 

with n — nn I a . Substituting these expressions into (7) 
yields the coupled system of linear algebraic expressions for 
the Hth Fourier component which can be partitioned and 
arranged in the following form: 



2. An actuating potential applied at the top surface, 



?ZT 



Va[nuJ : 



F -X s 



(9) 



III. Results 



Results are presented for simply supported hybrid beams. 
The thickness of piezoelectric layer is varied from 0.05/i to 

0. 30/z, where h is the total thickness of the beam. The beam 
configuration is as shown in Fig. 2. Beam (a) is a symmetric 4- 
ply laminate of equal thicknesses of material 1 with orientation 
[079079070°] while beam (b) is a sandwich beam with graphite- 
epoxy (material 2) faces and a soft core with thickness 8/z ; , 
where, /^is thickness of face. Beam (c) is an anti-symmetric 4- 
ply laminate made up of material 3 with lay-up [907079070°]. 
The total thickness of the beam is constant. The material 
properties are given in Table 1 . 

The results are plotted for uniformly distributed load as 
well as sinusoidal loading conditions in Figs. 3 through 5. 
Two load cases are considered: 

1. Uniform pressure p 2 = -p on top. 

2. Uniform applied potential f- = ^on the top. 
Similarly for sinusoidal loading conditions: 

1. A sinusoidal pressure on top surface, 

P \ =p sinixx/a)- 



Figure 2. Beam Configurations 

The results are non-dimensionalised with S = a/h, d T = 
374 x 10 12 CN Y T = 10.3 GPa for beam (a) and (c) and 6.9 GPa 
for beam (b): 

1. w = mwY T lhS i p ' CTx=C7 ' ISlp ^ 

2. w = \0wl S 2 d T <j> Q ,a x =(jh/WY T d T <j> , 

D : =Dh/l00Y T d;i- 



Table I. Material Properties 



Property 


Units 


Material 1 


Material 2 


Material 3 


Core 


PZT-5A 


Y, 


GPa 


181 


172.5 


181 


0.276 


61 


Y 2 


GPa 


10.3 


6.9 


10.3 


0.276 


61 


Y, 


GPa 


10.3 


6.9 


10.3 


3.45 


53.2 


G 12 


GPa 


7.17 


3.45 


7.17 


0.1104 


22.6 




GPa 


2.87 


1.38 


2.87 


0.414 


21.1 


G 31 


GPa 


7.17 


3.45 


7.17 


0.414 


21.1 


v ]a 




0.25 


0.25 


0.28 


0.25 


0.35 






0.25 


0.25 


0.28 


0.02 


0.3 8 


v „ 




0.33 


0.25 


0.33 


0.02 


0.38 


dm = dm 


pm/V 
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pm/V 
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dis = du 


pm/V 










584 


<1U = 122 


nF/m 










15.3 
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The effect of the beam lay-up is also studied by obtaining 
the results for symmetric [074574570°] as well as anti- 
symmetric angle ply beams [457-457457-45°]. The results 
are also obtained for a beam with piezo-layer attached to 
both top and bottom surfaces. These results are plotted in 
Figs. 6 through 9. 

Conclusions 

An efficient coupled electromechanical model has been 
presented for static analysis of hybrid beams with embedded 
or surface bonded piezoelectric layers under general 
electromechanical loading. This model combines a third order 
in-plane displacement and a layerwise linear electric field such 
that continuity of transverse shear stress at the interface of 
the layers is preserved with shear free conditions at top and 
bottom surfaces. The no. of primary variables is same as in 
FSDT and hence the computational cost is independent of 
no. of layers in the laminate. 

With increase in thickness of piezoelectric layer, the 
deflection of the beam increases and the normal stress at the 
top surface decreases as shown in the figures. It also reveals 
that the piezoelectric coupling effect increases with increase 
in piezo-layer thickness to overall thickness ratio. It is also 
observed that the deflection for anti-symmetric angle-ply 
laminates is more than that of cross ply and symmetric angle- 
ply laminates. The stresses increase for load case 1 while 
they decrease for load case 2. For beams with piezo-layers on 
both surfaces, deflection increases as the thickness of the 
piezo-layers is increased. 

Overall, the presented theory is well suited to model 
sensory as well as active response of smart composite beams 
for all applications. 
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Figure 3. Variation of }\> , (J and J) in beam (a) with thickness 
of piezo-layer for uniformly distributed loading 
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Figure 4. Variation of W , (J and £).in beam (c) with thickness of piezo-layer for sinusoidal loading 
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Figure 6. Variation of y\? and (7 in beam (a) with thickness of piezo-layer for uniformly distributed load case 1 
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Figure 7. Variation of W and (7 in beam (a) with thickness of piezo-layer for uniformly distributed load case 2 
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Figure 8. Variation of W , (J and £).in beam (a) with thickness of piezo-layer for sinusoidal loading 
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Figure 9. Variation of W , (J and £).in beam with thickness of piezo-layer on both surfaces for sinusoidal loading 
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